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Abstract 



We present an ab initio pseudopotential local density functional calcu- 
lation for stoichiometric high-T c cuprate YBa2Cu3C>7 using the plane-wave 
basis set. We have overcome well-known difficulties in applying pseudopo- 
tential methods to first-row elements, transition metals, and rare-earth ma- 
terials by carefully generating norm-conserving pseudopotentials with excel- 
lent transferability and employing an extremely efficient iterative diagonal- 
ization scheme optimized for our purpose. The self-consistent band struc- 
tures, the total and site-projected densities of states, the partial charges and 
their symmetry-decompositions, and some characteristic charge densities near 
Ep are presented. We compare our results with various existing (F)LAPW 
and (F)LMTO calculations and establish that the ab initio pseudopotential 
method is competitive with other methods in studying the electronic structure 
of such complicated materials as high- T c cuprates. 

PACS numbers: 74.72.Bk, 74.25.Jb, 71.25.-s, 71.20.-b 
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I. INTRODUCTION 



The discovery of superconductivity above 30K in La2_x.Baa.CuO4 by Bednorz and 
Miiller |l] has stimulated intensive studies on various cuprates as well as the mechanism 
underlying the phenomenon of high-T c superconductivity A wide variety of speculations 
about its microscopic origin have been raised with the sole consensus of the apparently im- 
portant role of the Cu02 layer (and, of course, the pairing of charge carriers) which is the 
basic building block of high- T c cuprates. Even for the normal (above T c ) state of the high- 
ly cuprates, the controversy as to whether it is a Fermi liquid or some other novel states 
has not been settled. The reason for the complication is that many of the cuprate supercon- 
ductors have sister compounds which are strongly correlated insulating antiferromagnets. 
Among novel points of view explaining anomaly in normal states of high- T c cuprates are 
the Luttinger liquid theory, the marginal Fermi liquid theory, the van Hove scenario, the 
spin fluctuation theory, the extended Hubbard model, the t- J model, and the anyon theory, 
to mention a few. On the other side of such efforts, there has been a more conventional 
approach dealing with the full Hamiltonian and the realistic atomic structure in the belief 
of the Fermi liquid theory; the first principles density functional theory |J with the local 
density approximation (LDA) || for the exchange- correlation energy is such an example. In 
spite of its shortcomings in the magnetic insulating phase, the various state-of-the-art elec- 
tronic structure calculations have proven to be successful in predicting from first principles 
a number of properties of metallic cuprates such as the lattice constant, atomic positions 
within the unit cell, phonon frequencies, structural instabilities, electric field gradients, and 
most importantly, the detailed energy band structure and the Fermi surface (FS) Q. Such 
calculations show the existence of a half-filled pcia-antibonding band and the resulting cylin- 
drical FS per Cu0 2 plane as the distinctive electronic feature of the high- T c cuprates || . 

Most of the existing electronic structures for high- T c cuprates have been calculated with 
use of the all-electron localized basis methods (LMTO, FLMTO, LAPW, or FLAPW) or 
the pseudopotential mixed basis method. Here, we want to emphasize the merits of the 
plane- wave basis set, which above-mentioned methods do not share. The mathematical 
formulation and the implementation of the numerical codes are particularly simple with 
the plane- wave basis set. The accuracy of calculation is easily controlled by varying the 
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kinetic energy cutoff. Plane-waves form a basis set independent of the ionic positions, 
allowing for an unbiased uniform description of the system and a much simpler calculation 
of quantum mechanical (Hellman-Feynman) forces. Furthermore, plane waves enable us to 
use the fast Fourier transform so that the iterative diagonalization methods IHTOll or the 
ab initio molecular dynamics JTT]^T~3[] can be applied efficiently In conjunction with the 
pseudopotential formulation [f[4l , fL5[| where the rapid oscillations of valence wavefunctions in 
the core region are systematically smoothed out with their norms kept conserved, the plane- 
wave basis sets are quite successful in the calculation of the electronic structure of most 
semiconductors and some metals. However, until recently it was believed that materials 
including first-row elements, transition-metal atoms, and rare-earth elements are not suitable 
for the application of the pseudopotential plane-wave method. With the development of 
efficient methods for diagonalizing large Hamiltonian matrices represented in a plane-wave 
basis set |6|-13j, together with relatively smooth pseudopotentials p6|Ji7[1 , it is now possible 
to perform the electronic structure calculations for systems containing such atoms. 

In this paper, we present an ab initio pseudopotential density-functional calculation of 
YBa 2 Cu 3 7 using the plane-wave basis set. YBa 2 Cu 3 7 _5 (0 < 5 < 1) is perhaps the 
most extensively studied compound among the cuprate superconductors. A unique feature 
of this compound is that the three copper atoms in the primitive cell play two different 
roles. Two copper atoms (Cu2) participate in two conduction layers of CuC>2 separated by 
yttrium atoms. The third copper atom (Cul) has an unusual coordination with oxygen 
atoms, forming a one-dimensional CuO chain structure along the b direction. The stoi- 
chiometric compound YP^CuaOy is particularly interesting from the standpoint of band 
theorists. It crystallizes in the relatively simple structure (with 13 atoms in the orthorhom- 
bic primitive cell) allowing for the computation with a less severe load compared with Bi- 
or Tl-based cuprate superconductors. More importantly, it is a metal instead of an anti- 
ferromagnetic insulator in its normal state in contrast to other stoichiometric cuprates like 
La2Cu04. Aside from the problem in handling the strong correlation among electrons, the 
local density functional (LDF) calculation is applicable to this system. Furthermore, clean 
surfaces are available in this compound so that high resolution spectroscopic experiments 
can be performed and the theoretical results can be compared directly with experiments. 
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The quasiparticle band near the Fermi level and the FS topology from the ARPES [18 



ACAR [|19|| , and dHvA [20| experiments are in good accordance with those obtained from 



the existing LDF calculations [pTH 29[|. This fact is sometimes regarded as supporting evi- 
dence for the Fermi liquid picture of the normal state of high- T c cuprates though the issue 
is still controversial. Since the main purpose of the present work is to show the feasibil- 
ity of the ab initio pseudopotential plane- wave method for studying high-T c materials in 
comparison with other methods and since the present paper is the first such report to our 
knowledge, we will exhibit below a rather extensive amount of figures obtained through the 
fully-converged pseudopotential calculation, but avoid promoting a particular theory for the 
high- T c superconducting mechanism. 

II. COMPUTATIONAL METHOD 

We generate soft norm- conserving pseudopotentials using the scheme of Troullier and 



Martins JT6[, which gives a very fast convergence of the total energy with respect to the 
basis-set size. (It is possible to use ultrasoft pseudopotentials which require a much 
smaller basis set size at the price of abandoning the norm- conserving property, but we do 
not adopt that scheme here.) To enhance the transferability of the pseudopotential and 
guarantee correct scattering properties of the pseudopotential, we include some semi-core 
orbitals into the valence shell and use ionized configurations for the atoms Y and Ba. The 
configuration and the cutoff radii r^'s with which the pseudopotentials are generated are 
listed in Table |I[ All the pseudopotentials, except for the oxygen pseudopotential, are gen- 



erated semi-relativistically. The partial core correction scheme |3l| is employed to overcome 
the problem associated with the nonlinearity of the exchange-correlation functional. The 
Ceperly-Alder correlation [3^] is used in the parametrized version of Perdew and Zunger 



These pseudopotentials are then cast into the fully-nonlocal separable form of Kleinman and 
Bylander [Q with use of s-locality to avoid the ghost states [|3^,[!(J. The norm- conserving 
requirement results in the matched first energy derivative of the logarithmic derivative of 
the wavefunction, logi?/(r)| e=eijT . =r . cP between the pseudo radial wavefunction and the 
all-electron one, where Ri(r) and ei are the radial wavefunction and the valence orbital en- 
ergy for angular momentum I. The pseudopotentials optimized by us with the above scheme 



show scattering properties acceptable in a much wider energy range than enforced by the 
norm-conserving condition. 

We then perform the plane-wave band structure calculations within the LDA. The 



single-particle Kohn-Sham equation is solved by diagonalizing the Hamiltonian matrix |37 
Diagonalization is achieved by using the block-Davidson method |5],|7| with the modified 
Jacobi relaxation |8|-|T0lj. At this point, we want to mention the superiority of the block- 



Davidson method with the modified Jacobi relaxation to two other competing methods, 
the DIIS- ritzit |^-|ID| and the preconditioned conjugate gradient [0,0 methods, in dealing 
with systems with relatively localized states. We have tested these three diagonalization 
schemes by studying the cubic ZnS with variable submatrix sizes. For relatively small 
sub-matrix sizes, the block-Davidson method still gives correct results, whereas the latter two 
methods frequently fail to do it. The discrepancy is due to the fact that the block-Davidson 
method updates the eigenvectors simultaneously while the latter two adopt a band-by-band 
updating scheme. Actually, there exists a crossover behavior in relative merits among the 
three schemes. If the number of desired states is very large, the computational load in 
diagonalizing the Hamiltonian projected to the expansion space is too severe with the first 
method and the latter two methods with a sufficiently large submatrix is preferable to the 
first one. As the number decreases, however, the computational effort of the first method 
in the diagonalization of the expansion- space-projected Hamiltonian is greatly reduced and 
the block-Davidson method is superior to the other two. 

III. RESULTS AND DISCUSSION 

The calculation for the stoichiometric YE^CuaOy are performed with use of the recently 
obtained structural parameters at 10K p8||; a = 3.8196A, b = 3.8813A, c = 11.64028A, z Ba , 



= 0.18388, z Cu2 = 0.35457, z Q2 = 0.37823, z Q3 = 0.37813, and z Q4 = 0.15976. The atomic 
positions in lattice unit are Y(|, |, |), Ba(~, |, ±z Ba ), Cul(0, 0, 0), Cu2(0, 0, ±z Cu2 ), Ol(0, 
§, 0), 02(|, 0, ±z 02 ), O3(0, §, ±z Q3 ), and O4(0, 0, ±z ot ). The cutoff energy of 72.25 Ry 
and 147 k-points in the irreducible Brillouin zone (IBZ, 1/8 of the Brillouin zone) with the 
linear tetrahedron method p9| , f4"0f are used to obtain the self-consistency. The self-consistent 
energy bands are plotted in Fig. [TJ along the symmetry directions in the IBZ. The same band 



structure is shown in expanded scale in Fig. |2| for the energy range from —1.0 eV to 3.0 eV. 
There are four partially filled bands crossing Ep. Two nearly half-filled bands are derived 
from the two CuC>2 planes and have the usual two-dimensional character of Cu2(3<i r 2_, / 2)- 
02(2p x )-03(2p y ) er-antibonding. The other two are derived from the CuO chain and show 
negligible dispersions along S-Y and R-T in contrast to appreciable dispersions along X-S and 
U-R. The two CuO chain bands show completely different bonding characters. The nearly- 
empty dispersive band has the character of Cul(3<i y 2_ 2 2)-01(2p 2/ )-04(2p 2 ) a-antibonding 
and the almost filled flat band has the Cul(3d yz )-01(2p z )-04:(2py) 7r-antibonding character. 
Four ^-crossing bands form a complicated FS's shown in Fig. one stick, two cylinders 
along the c direction, and one sheet perpendicular to the T-Y direction. There are many 
published literatures presenting either the energy band structure or the FS's of YBa2Cus07, 
all of which resort to localized basis methods. Among them, the first report was made by 



Yu et al. pl] , |22]| and illustrates an overall feature of the band structure of YB^CusOy, but 
shows a somewhat unconverged result — the existence of a hole-like cylindrical FS along 



Y-T. Later Massidda |24] assigned the incorrectness to the small number (10) of k-points 
and corrected it with use of a larger number (32) of k-points. The second report by Krakauer 



et al. PBj shows more converged results than Yu et al., but the interaction between plane 



bands and the chain band are incorrectly described as will be discussed below. Results 



of the subsequent calculation by the same group (Pickett et al. [23|) are different from its 
preceding work; the band-band interaction is correctly described but the chain derived band 
along T-Y, which is slightly unfilled in Krakauer et al. |[23|| , lies below Ep. Andersen et a/.'s 
FLMTO results [p7| are indistinguishable from Pickett et a/.'s. There is a calculation done 



by using the pseudopotential mixed-basis method |28|| , but the resulting energy bands show 
poor convergence. Judging from the band structure and the calculated FS's, the present 
pseudopotential plane-wave calculation shows the most converged results with accuracy 
comparable with Pickett et al. and Andersen et al. ||27|| . (Though the FS topology is 
different for the large plane-derived band, it is fragile with respect to the variation of Ep 
on a scale as small as ~ 0.01 eV.) The present results are also in good agreement with 
experiments [|i~8|-|2"0| . 

The bonding nature of the four partially filled bands can be clarified by examining their 
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charge density contour plots at some selected k-points. The charge densities of the states 
at k = (|, 7, 0), a half way from T to S, for the lower and the upper Cu0 2 bands are shown 
in Figs. |](a) and (b) respectively. Both of them reflect the Cu2(d x 2 _ y 2)-02(p x )-03(p y ) a- 
antibonding hybridization characteristic of the high-T c cuprates. In detail, however, their 
symmetries are different. The wavefunction of the lower CuC>2 band is the symmetric linear 
combination of the two CuC>2 plane states residing above and below the Y atom plane, and 
the upper CuC>2 band is the antisymmetric one, which can be deduced from the comparison 
of the smallest- value contour lines in Figs. |](a) and (b). Furthermore, their parity with 
respect to the chain plane reflection is even for the lower band and odd for the upper one, 
respectively. Since the wavefunctions of the nearly-empty chain band are even with respect 
to the reflection in the plane of the chains perpendicular to c axis, it follows that for k z = 
the chain band does cross the odd (upper) plane band while hybridizing with the even (lower) 
plane band. For k z = -, the situation is reversed, i.e., the nearly-empty chain band interacts 
with the upper plane band but does not with the lower one. Such symmetry properties and 
the interaction between the plane bands and the nearly-empty chain band are mentioned in 



Pickett et al. |25| and fully described in Mazin et al. 

The charge densities of the nearly-empty chain band are shown in Fig. [5| at k = (~, |, 0) 
and (4, 4, 0). The Cul(c? y 2„ 2 2)-01(pj / )-04(^) a-antibonding character and the local four- 
fold coordination of the Cul atom are clearly identified from them. The absence of the band 
dispersion along S-Y or R-T can be ascribed to the strong one-dimensional character along 
the b axis. As the energy rises towards the S point the antibonding character of the Cul-Ol 
bond increases while that of Cul-04 is almost unaffected. The shorter Cul-04 bond length 
of 1.8596A compared with the Cu2-04 interatomic distance of 2.2676A is consistent with 
the bonding characters identified above. The Cul-Ol-04 pcia-antibonding band is almost 
empty, while its bonding counterpart is completely filled. It is therefore likely for 04 atoms 
to be tightly bonded to Cul atoms. On the other hand, the Cu2 atoms have no ligand O 
atom on the Y plane side. The Cu2(d 3z 2_ r 2) orbitals are only weakly bonded to 04 and 
both the bonding and antibonding bands are almost filled, resulting in a weaker Cu2-04 
bond and a longer Cu2-04 interatomic distance. 

The contour plots of the charge density for the fourth (almost filled) band are shown in 
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Fig. H at k = (~, ~, 0) and (|, |, 0). A remarkable feature is that the main contribution 
to this band comes from 01(p z ) and OA(p y ) orbitals and relatively little charge resides at 
the Cu, 02, and 03 atoms. The character at S is 01(p z )-04:(py)-Cul(dy Z )-Cu2(dy Z ) ir- 
antibonding and the charge is depleted on the zx-plane. The states near k =(|, 7, 0) have 
the 01(p z ), 04(p x , p y , p z ), Cul(d yz , d zx ), Cv2(d z 2_ x 2), and 02(p z ) orbital characters. Since 
these chain-related states lie in the narrow energy region containing Ep, they are expected 
to be strongly affected as the oxygen content on the CuO chain is reduced and the hole 
doping concentration decreases, implying that T c should be influenced significantly by these 
states. 

The muffin-tin radii of 2.80, 3.20, 1.95, and 1.65 a.u. for Y, Ba, Cu, and O atoms re- 
spectively are used in calculating the partial charges and the partial (or site-projected and 
symmetry-decomposed) densities of states (PDOS). The total density of states (TDOS) 
at Ep is 1.41 states/(eV Cu). For comparison, Massidda et al. |21| obtained N(Ep) = 
1.13 states/(eV Cu) and Krakauer et al. [^3J 1.38. The PDOS plots in Fig. |7| show that 
the Cul 01 04 sates prevail in the energy range from 0.5 eV below Ep to Ep. The PDOS 
in an expanded scale in Fig. |] indicates a van Hove singularity at ~ 0.03 eV below Ep 
with a clear Cul-Ol-04 hybridization. The importance of the van Hove singularity to 
the high-T c mechanism has been strongly advocated in some theoretical studies [[|l]]. The 
partial charges and their symmetry decompositions are listed in Table [H| (The partial 
charges and their symmetry decompositions are presented and fully analyzed in Schwarz 
et al. [p6| .) The relatively small amount of 01(p y ), 04(p z ), Cul(d x 2_ y 2), and Cul(G? 3z 2_ r 2) 
partial charges are due to the nearly-empty dispersive band of Cul(c? 2/ 2_ z 2)-01(pj / )-04(p z ) 
hybridization. In fact, with use of a different choice of symmetry decomposition includ- 
ing d y 2_ z 2 symmetry, the partial charge of Cul(<i y 2_ z 2) is further reduced to 1.30, which 
is smaller than the partial charge of Cul^^.^), 1.38. The relatively small amount of 
02(p x ), 03(p y ), and Cu2(d a .2_ J/ 2) partial charges compared with the other 02(p), 03(p), 
and Cu2(c/) partial charges is attributed to the characteristic pda- antibonding hybridization 
[Cu2(<i r 2_, y 2)-02(p :E )-03(]9j / )] of the superconducting cuprates. The smaller partial charge 
of Cul (3^2.22) than that of Cu2(3c? :r 2„ y 2) (1.30 versus 1.42) is in accordance with the ex- 
istence of the nearly-empty chain band and two nearly half-filled plane bands. The partial 
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charges of the Cu(4p) states are relatively small, but their magnitude can serve as a measure 
of the bond strength between a copper atom and its neighboring atoms. The Cu(4p) partial 
charges do not originate from on-site Cu(4p) states, which lie at much higher energies, but 
represent the intrusion of wavefunctions centered at neighboring (off-site) atoms into the 
copper muffin-tin sphere. They originate mainly from the surrounding 0(2p) wavefunctions 
and their magnitude depends on the distance from copper to the neighboring oxygen sites. 
Since the Cul-04 bond length is shorter than that of Cul-Ol, the partial charge of Cul(jo 2 ) 
is larger than that of Cul(p y ). The Cul(p x ) partial charge is the smallest due to the absence 
of the oxygen neighbors in the a direction. For the plane copper atoms, the longer distance 
to 04 than to 02 or 03 manifests itself in the smaller Cu2(p z ) partial charge than that of 
Cu2(p x ) or Cu2(p„). 

IV. SUMMARY 

We have applied, for the first time, the ab initio pseudopotential plane-wave method 
to the electronic structure calculation of the high-T c cuprate YE^CusOy within the LDF 
formalism. Soft pseudopotentials of optimized transferability are generated through the 
Troullier-Martins scheme and cast into the Kleinman-Bylander type fully-separable form. 
The block-Davidson algorithm with the modified Jacobi relaxation operator is then used to 
solve the single-particle Kohn-Sham equation. The energy bands, the Fermi surfaces, the 
charge densities, the total DOS, the partial DOS, and the symmetry-decomposed partial 
charges obtained here are in good agreement with both the existing theoretical calculations 
using (F)LMTO and (F)LAPW methods and the available experimental facts. This success- 
ful application to YBa 2 Cu 3 7 has demonstrated that one can deal with systems containing 
first row elements, transition metals, or rare-earth elements with acceptable computational 
efforts via the ab initio pseudopotential plane-wave method. 
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TABLES 

TABLE I. The reference configurations and core radii. For Y and Ba atoms, +3 and +2 ionized 
configurations are used respectively to get optimal pseudopotentials. 



atom 


Y 


Ba 


Cu 


O 


configuration 


4 s 2 4p 6 4d i 


5s 2 5p 6 5d° 


4s Hp 3d 10 


2s 2 2p 4 


^cs 


1.19 


1.65 


2.06 


1.15 




1.35 


2.00 


2.30 


1.45 


r c d 


2.49 


2.48 


2.06 





TABLE II. The symmetry-decomposed partial charges (in electrons) of the valence states of 
YBa 2 Cu 3 7 . 



atom 


s 


X 


y 


z 


P 


xy 


yz 


zx 


x 2 - y 2 


3z 2 - r 2 


d 


total 


Y 


0.086 


0.069 


0.066 


0.070 


0.204 


0.048 


0.136 


0.146 


0.137 


0.052 


0.521 


0.811 


Ba 


0.062 


0.135 


0.123 


0.120 


0.377 


0.092 


0.039 


0.058 


0.037 


0.060 


0.286 


0.725 


Cul 


0.194 


0.028 


0.059 


0.079 


0.167 


1.836 


1.856 


1.837 


1.644 


1.379 


8.552 


8.912 


Cu2 


0.175 
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Figure Captions 

FIG. 1. The energy band structure of YBa2Cus07. 

FIG. 2. The energy band structure of YBa 2 Cu 3 7 on (a) k z = plane and (b) k z = 
- plane in expanded scale. The inset in (b) is the band structure in much more 
expanded scale illustrating the band crossing between the chain band and the 
lower plane band as well as the band repulsion between the chain band and the 
upper plane band. 

FIG. 3. The FS's of YBa2Cus07 derived from the four partially filled bands. The solid 
line is the stick-like FS along S-R from the almost filled chain band, the short- 
dashed line indicates the one- dimensional FS arising from the nearly-empty chain 
band, the long-dashed line is the cylindrical FS from the upper Cu02 band, and 
the medium-dashed line from the lower Cu02 band. 

FIG. 4. The charge density contours of YBa2Cus07 for states at k = (|, |, 0) of (a) 
the lower Cu2-02-03 plane band and (b) the upper Cu2-02-03 plane band. 
Contour values run from 0.0002 to 0.0102 in steps of 0.002 e/unit-cell. The 
symbols # and ♦ represent Cu and O atoms respectively. 

FIG. 5. The charge density contours of YBa2Cu30y for the nearly-empty Cul-Ol-04 
cr-antibonding chain band at (a) k = (|, |, 0) and (b) k = (|, |, 0). Contour 
values run from 0.001 to 0.011 in steps of 0.002 e/unit-cell. 

FIG. 6. The charge density plots of the almost filled chain-related band of YBa2Cus07 
at (a) k = (|, 0) and (b) k = (i, |, 0). Contour values are the same as in 
Fig. |. 

FIG. 7. The total and the partial densities of states of YBa2Cu307. (a) The TDOS, 
(b) the PDOS of chain-related atoms, and (c) the PDOS of plane-related atoms. 

FIG. 8. (a) The TDOS, (b) the PDOS of chain-related atoms, and (c) the PDOS of 
plane-related atoms of YBa2Cus07 near Ep in expanded scale. 
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